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RESONANCE  BETWEEN  BETATRON  AND  SYNCHROTRON 
OSCILLATIONS  IN  A  FREE  ELECTRON  LASER: 

A  3-D  NUMERICAL  STUDY 


INTRODUCTION 

In  this  paper,  we  will  examine  the  resonance  between  the  synchrotron  wavenumber, 
i.e.,  the  wavenumber  associated  with  the  bounce  of  the  electron  in  the  ponderomotive 
potential  well,  and  twice  the  betatron  wavenumber  [1-5].  Betatron  oscillations  are  caused 
by  the  transverse  gradient  of  the  wiggler  [6-12].  There  are  a  number  of  different  sources 
for  this  resonance.  In  this  paper,  we  will  concentrate  on  the  one  caused  by  the  Gaussian 
resonator  radiation  field  first  discussed  by  M.  Rosenbluth  in  Ref.  [l].  Another  source  of 
the  resonance,  described  in  detail  in  Ref  [5],  requires  a  tapered  wiggler,  but  no;  a  3-D 
radiation  field. 

The  wavefronts  of  the  Gaussian  TEMyo  resonator  radiation  are  spherical;  the  curva¬ 
ture  Ls  a  function  of  radial  and  axial  position.  As  electrons  execute  betatron  oscillations  in 
the  transverse  direction,  the  electrons  are  forced  to  sample  varying  wavefronts  or  a  varying 
phase  of  the  radiation  field.  The  effect  of  that  on  the  axial  electron  motion  is  equivalent 
to  varying  the  radiation  frequency,  which  results  in  a  forced  driving  term  in  the  modified 
pendulum  equation  with  a  wavenumber  of  2kp,  where  k$  is  the  betatron  wavenumber. 

The  amplitude  of  the  TEMU<,  resonator  radiation  field  has  a  Gaussian  profile  in  the 
radial  direction.  As  electrons  execute  betatron  oscillations  in  the  transverse  direction, 
the  electrons  are  forced  to  sample  varying  radiation  amplitudes.  In  an  FEL  employing 
efficiency  enhancement  schemes,  we  find  that  the  amplitude  variation  also  contributes 
forced  driving  terms  in  the  modified  pendulum  equation  with  wavenumber  2 kj. 

Betatron  oscillations,  in  a  wiggler  whose  magnetic  field  amplitude  Bv.  is  tapered, 
result  in  a  spatially  dependent  axial  electron  momentum.  The  wavenumber  of  the  axial 
momentum  variation  is  also  2 k?  (see  Ref.  [5]). 

When  2 kg  becomes  approximately  equal  to  the  synchrotron  wavenumber  Ks,  tl  e  phase 
of  the  electrons  in  the  ponderomotive  potential  well  can  become  unstable  and  electrons  can 
become  detrapped  from  the  buckets.  Since  the  principle  of  efficiency  enhancement  schemes 
[13-1G]  for  the  FEL  is  based  on  trapping  the  electrons,  the  loss  of  trapped  electrons  will 
Manuscript  approved  October  17,  1984. 
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reduce  the  gain.  In  an  FEL  oscillator,  a  slightly  reduced  gain  will  increase  the  start-up 
time  necessary  to  reach  saturation.  A  significantly  reduced  gain  per  pass  will  also  reduce 
the  final  saturated  field  amplitude  as  well, 
amplitude  as  well. 

In  this  paper,  we  will  consider  a  realistic,  arbitrarily  tapered,  linearly  polarized,  mag¬ 
netic  wiggler.  The  vector  potential  of  the  wiggler  is  expressed  as 

■Aj u  =  Aw [z}cosh(kw  (-)y)  cos(  /  kw{z  )dz  Je^,  (1) 

Jo 

where  Aw(z)  and  kw(z)  are  the  slowly  varying  amplitude  and  wavenumber.  We  will  con¬ 
sider  an  electron  beam  that  is  cold  but  possesses  a  finite  emittauce.  The  radiation  field 
initially  is  taken  to  be  a  Gaussian  T EAf0o  resonator  mode  with  slowly  varying  amplitude 
and  phase.  Modeling  of  an  exact  3-D  radiation  with  higher  order  modes  can  be  acheived 


with  minor  modifications. 


For  a  linearly  polarized  wiggler,  there  exists  a  constant  of  motion  in  the  x-direction, 


Px  Ax  —  Pxoi 

c 


where  the  constant  pxo  is  the  canonical  momentum.  To  obtain  the  betatron  orbits,  we  will 
take  Ax  a*  Aw  ■  ex,  and  px  c*  ^Ax  +  pxo-  The  particle  motion  in  the  x-direction  can  be 
obtained  by  integrating  the  momentum  in  the  x-direction, 


-  .  Po±{z) 

X  X  +  — — J-T  COS 

kw[z) 


Hkw(z)y)sin(J  kw(z')dz'), 


where  x  =  x0  +  0xoz,  fixo  =  Pxo/irn0c,  y  is  a  function  of  (y0,Pyo,z),  and  y0  are  the 
initial  transverse  coordinates,  and  pxo  and  pyo  are  the  initial  transverse  momentuins,  and 
0o±  =  |e|Au,/(7m0c2). 

The  particle  motion  in  the  y-direction  is  due  to  the  finite  z-component  of  the  wiggler 


We  will  assume  that  the  fast  oscillatory  terms  are  unimportant.  Replacing  vz  by  c  in  the 
appropriate  places,  we  find  that 


+  ^0a±(z)kw(z)  sinh(2 kw{z)y)  =  0.  (5) 

Taking  kw(z)y  <<  1,  this  equation  can  be  integrated  to  give 

V-Vp  (j^j)  kfi(z')dz'  +  (6) 

where  k #  =  fl0±kw/\/2,  Vp  =  {fyo  +  vlf1*  .  Pyo  =  Pyolkp[Q)vn0(? , 
and  <J>0  =  cos ~l{y0/vp)- 


RADIATION  FIELD 

The  radiation  field  can  be  written  in  terms  of  Gaussian  resonator  modes.  For  the 
purpose  of  illustrating  this  resonance,  only  the  TEMqq  mode  will  be  considered,  i.e., 

z,t)  c*  -Aoo(*)Goo(*,  y,  z)  exp(i(fcz  -  +c.c.,  (7) 


where  Aoo(-)  =  |A0o(^)|  exp(i'0Oo(^))  is  the  slowly  varying  self-consistent  complex  ampli¬ 
tude, 


Goo(*,y,  z)  =  goo  exp(«0Oo), 

,  x  wo  M*3  +  yan 

goo(*,V,z)  =  — j-rexp  I  - jj-r —  1  , 

w[z)  \  wl[z)  ) 

(X2  +y2  \ 

w2(zf) 


(8) 


f  =  (z  -  Lc)/z0  is  the  normalized  axial  distance,  w(z)  =  u>0(l  +  f2)1/2,  w  and  k  =  w/c 
are  the  frequency  and  wavenumber  of  the  radiation,  w0  is  the  waist  of  the  resonator  mode, 
z„  =  w*kl 2  is  the  Rayleigh  length,  and  Lc  is  the  axial  location  of  the  minimum  waist. 

Gaussian  decomposition  of  the  radiation  field  has  been  outlined  in  several  different 
contexts  (2, 17].  This  formulation  differs  from  Ref.  [17]  in  that  we  include  betatrcn  oscil- 


The  slowly  varying  complex  amplitude  of  the  radiation  held  is  governed  by 


da°o(z )  =  iaFa|u(3)e««oo(^)  ^cosh(fcr(;)tf)goo(i,y,j)c_,y^  ^ 


(9) 


where 


a  = 


vzo 


2k  c1  c 

'P{*o,Uo,Pio,Pyo ,4>o,z)  =  f  {k  +  ku,[z')  -  u)/vz)dz'  +0Oo{i,y,z)  +  <Poo(z)  +  tpo 

J  o 


is  the  phase  of  the  electron  in  the  ponderomotive  potential  well,  t[’0  is  the  initial  phase  at 

2  =  0, 


(('"))“  j  J  dx°  f*V°  J  dP10  J  dPyo {-)W{z0,y0,pXo,Pyo,  4’o ) 


is  the  average  over  all  electrons  in  the  transverse  direction  and  one  period  of  the  pon¬ 
deromotive  potential  wave,  W  is  the  initial  electron  distribution  function  such  that 
<  (l)  >=  1,  Ub  =  (4jr|e|2n0/m0)U2  is  the  plasma  frequency,  n0  is  the  peak  electron 
density,  F  =  <rb/irw2  is  the  1-dimensional  filling  factor,  <rb  is  the  area  of  the  electron  beam, 
aw  =  \e\Awjm0c2  and  a00  =  \e\A^o(m0c2 . 


EQUATION  OF  PHASE 


The  equation  of  phase  for  the  electron  entering  the  wiggler  with  the  initial  condition 
(x0,  Vo,Pxo,Pyo,  V’o)  in  the  ponderomotive  potential  well  can  be  written  in  the  form 


d2ip 


dz2 


=  -K20  coshf&w  ( z)y)e 


\  "  ,x>  /  Isintfc -si 


sin  0  —  sin  , 


(10) 


where 


K'o{z)  =  (^^***a*(z)|aoo(z)|) 


1/2 


is  the  synchrotron  wavenumber  of  the  electrons  traveling  exactly  along  the  z-axis.  The 
instantaneous  synchrotron  wavenumber  for  the  off-axis  electrons  is 

1/2 


K,{x0,y0,Pxo,Pyo,z)  =KS0(z )  (cosh{kw(z)y)exp(-^-^r-^ 


T(  \  -  _  k  da*  .  +  je| 

*  [  dz  472  dz  c  7s  m0c2  5s  J 

is  the  degree  of  taper  for  the  efficiency  enhancement  schemes,  which  can  be  achieved 
by  increasing  the  wiggler  wavenumber,  decreasing  the  amplitude  of  the  wiggler  vector 
potential,  or  by  applying  a  DC  accelerating  electric  field.  For  efficiency  enhancement,  the 
resonant  phase  has  to  be  positive  and  less  than  unity. 

In  the  betatron-synchrotron  instability  the  resonant  phase  is  forced  to  oscillate  at  the 
wavenumber  2kp.  If  the  wavenumber  of  the  electrons  about  the  separatrix  is  ^  2 kp)  then  an 
instability  can  occur.  As  electrons  execute  betatron  oscillations,  the  electrons  experience 
a  varying  radiation  amplitude  and  phase,  which  are  contained  in  the  second  and  the  third 
terms  on  the  right-hand  side  of  Eq.  (11).  The  tapering  of  the  wiggler  Bw  can  also  result  in 
an  oscillatory  driving  term,  i.e.,  the  fourth  term  on  the  right-hand  side  of  Eq.  (11).  This 
process  is  discussed  in  detail  in  Ref.  |5j.  If  the  oscillation  frequency  of  the  resonant  phase 
matches  the  frequency  of  the  electrons  going  around  the  resonant  phase,  the  synchrotron 
frequency,  then  the  phase  of  the  electrons  in  the  ponderomotive  potential  well  will  oscillate 
with  increasing  amplitude,  and  eventually  the  electrons  become  detrapped. 
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The  radiation  amplitude  variation  felt  by  the  electrons  not  only  results  in  a  spatially 
dependent  bucket  size,  but  also  in  forced  oscillations  of  the  resonant  phase,  the  second 
term  of  Eq.  (11),  where 


=  1  - 


l _  fi2  +  y2\ 

M~')y)expV  «’2U)  )' 


cosh(kw(:)y)  V  «’2(-)  ) 

Since  y 2  and  cosh(fcu,(s)y)  are  both  periodic  functions  of  $,  the  coefficient  of  the  second 
term  in  Eq.  (11)  is  also  a  periodic  function  of  $. 

The  variation  in  phase  felt  by  the  electrons,  due  to  the  betatron  oscillations,  con¬ 
tributes  the  third  term  on  the  right-hand  side  of  Eq.  (11).  Substituting  (6)  into  (8),  we 
obtain 


=  /„  +  fe  cos  *  +  /,  sin  *, 


fs{y0.PyoiZ)  =  ~2 


y I  M*)  [i-?2' 

U>2(*)  So  [l  +  f2. 


ft  \  1  ?2~3 

Jo[*o,Vo,Pxo,P*o n-J-  22(1  +  f2j2  2W*(Z) 


~2 ,  M»)  n 


3?  *■ 


w*(z) 


2&IOX  1  f2  ,  j2  . 

— rr?+M- 


We  note  that  each  term  on  the  right-hand  side  of  Eq.  (12)  is  highly  dependent  on  the 
initial  conditions  as  well  as  the  axial  position  z.  The  phase  variation  term  (l/K^cPOoo/ds2 
does  not  have  a  spatitial  variation  equal  to  2 kp,  because  the  coefficients  A's,  /s,  /c,  and  f0 
are  functions  of  z.  Hence,  the  effectiveness  of  the  phase  variation  term  in  driving  the 
betatron-synchrotron  instability  is  reduced. 

The  strongest  effect  of  the  resonance  on  the  gain  is  not  at  Kso  =  2^(0).  becfuse  the 
synchrotron  wavelength  associated  with  electrons  undergoing  betatron  oscillation  i>  longer 
than  that  associated  with  electrons  traveling  exactly  on  axis.  This  is  a  consequence  of 


the  larger  radiation  field  on  axis.  In  addition,  the  synchrotron  wavelength  of  the  electrons 
initially  close  to  the  resonant  phase  iI’r  is  2 ir/Ks.  The  synchrotron  wavelength  of  the 
electrons,  that  were  initially  trapped  further  away  from  the  resonant  phase,  is  longer  than 
2itjK9.  Hence,  the  more  accurate  resonant  condition  is 

6 >  KSo  >  2 k*.  (13) 

NUMERICAL  EXAMPLES 


We  present  numerical  results  illustrating  various  properties  of  this  instability.  The 
linearly  polarized  wiggler  amplitude  is  5^.(0)  =  3  kG,  and  the  wavelength  is  ^.(0)  =  2.73 
cm.  The  wiggler  length  is  Lw  =  150  cm,  the  minimum  waist  is  located  at  z  =  75  cm,  and 
the  electron  beam  energy  is  20.8 MeV.  The  electron  density  on  axis  is  n0  =  2.5xlOu  cm-3, 
which  corresponds  to  a  plasma  frequency  u of  2.8  x  1010  sec-1.  Radiation  wavelength  is 
10  /< m,  with  a  Rayleigh  length  zD  =  62.5  cm  and  spot  size  wD  =  0.14  cm.  The  betatron 
wavenumber  is  kp  =  0.03  cm-1  and  the  betatron  wavelength  is  Lp  =  2ir/kp  =  21C  cm. 

A  Gaussian  electron  distribution  function  is  chosen: 


W  —  C  exp 


where  C  is  the  normalization  constant,  and  pxo  =  0. 

Here  we  present  results  where  kw  is  tapered  and  aw  is  held  constant.  The  taper  of  kw 
is  linear,  i.e., 

kw(z)  =  &a/(0)(l  +  6z), 

where  6  =  1/LC,  and  Lc  is  the  length  the  wavenumber  increases  by  &u;(0)-  For  .i  linear 
variation  of  k sin  V’r0  is  not  a  constant.  The  constant  6  is  chosen  such  that  sin  i  r0[z  = 
Lu.)  =  0.4.  The  parameters  that  will  be  varied  are  the  betatron-synchrotron  waveiumber 
mismatch  ratio 

_  Kso{,z  ~  Lc) 

P~  »/>(0)  ’ 


and  the  radius  of  the  electron  beam  reb ■ 


Figures  16  show  the  results  obtained  with  p  =  2.0  and  reb  =  «'o-  For  these  parameters, 
the  second  term  on  the  right-hand  side  of  Eq.  (11)  is  the  largest  oscillatory  driving  term. 
For  electrons  with  the  same  initial  conditions  (x0,y0,py0)  but  -;r  <  <  tt.  Figs.  1,  3,  5 

are  plots  of  f  versus  dii'jdz  and  Figs.  2,  4,  6  are  plots  of  v  versus  z.  Each  curve  represents 
the  history  of  an  electron  in  the  phase  space  from  ;  =  0  to  c  =  Lu.. 

Figures  1  and  2  are  plots  of  phase  for  electrons  traveling  exactly  down  the  a>ris,  i.e., 
x u  =  0,  y0  =  0,  and  py0  =  0.  Electrons  that  are  initially  trapped  remain  trapped. 

Figures  3  and  4  are  plots  of  phase  for  electrons  with  initial  conditions  i0  —  w0,  y0  =  0, 
and  pyo  =  0.  These  electrons  do  not  execute  betatron  oscillations.  Since  they  are  on  the 
edge  of  the  radiation  beam,  the  bucket  size  is  reduced,  and  only  40^i  of  the  electrons  are 
initially  trapped,  these  trapped  electrons  remain  trapped  to  the  end  of  the  wiggler. 

Figures  5  and  6  are  phase  plots  for  electrons  with  initial  conditions  x0  —  0,  yc  =  w0. 
and  pyo  =  0.  If  these  electrons  did  not  undergo  betatron  oscillations,  their  phase  space 
diagrams  would  be  identical  to  those  in  Figs.  3  and  4.  We  notice  that  25!^  of  the  electrons 
become  detrapped  at  ;  ~  100  cm  jr/fc^O).  Only  20°t  of  the  electrons  remain  trapped 
at  the  end  of  the  wiggler. 

The  numerical  results  confirm  that  the  synchrotron  wavelengths  associated  with  Figs. 
3  6  are  much  longer  than  those  of  Figs.  1-2,  because  the  radiation  amplitude  felt  by 
the  electrons  on  axis  is  larger  than  the  the  average  radiation  field  felt  by  the  electrons 
undergoing  betatrou  oscillations. 

Figure  7  is  a  plot  of  normalized  relative  amplitude  gain, 

C  -  g/r«fc 
*  f/r*6|rrt- 0 

versus  reh/w0  for  p  =  1.4,  p  =  2.0,  and  p  =  2.9.  The  solid  curves  are  the  results  with 
betatron  oscillations,  and  the  dashed  curves  are  the  results  without  betatrou  oscillations. 
The  results  are  almost  identical  for  reb/vc  <  0.8.  There  is  a  small  reduction  of  the  gain 
when  reb/ir0  =  1  for  p  =  2.0  and  p  =  1.4.  The  reduction  is  small  because  the  distribution 
function  associated  with  the  electrons  that  take  part  in  the  betatron-synchrotron  instability 
is  small. 
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In  conclusion,  we  identified  three  different  sources  responsible  for  the  betatron-  syn¬ 
chrotron  instability:  the  radiation  phase  front  curvature,  the  transverse  radiation  ampli¬ 
tude  variation  in  an  FEL  with  efficiency  enhancement  schemes,  and  the  tapering  of  the 
magnetic  wiggler  field  amplitude.  For  the  parameters  considered  in  this  paper,  the  radi¬ 
ation  amplitude  variation  is  found  to  be  the  largest  driving  term.  For  short  wigglers  of 
this  example,  i.e.,  Lw  ~  1.5 and  reb/w0  5:  0-8,  the  effect  of  betatron-synchrotron 
instability  on  the  gain  is  negligible. 
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